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Two-component fractional quantum Hall systems are providing a major motivation for a large 
section of the physics community. Here we study two-component fractional quantum Hall systems 
in spin-polarized half-filled lowest Landau level (filling factor 1/2) and second Landau level (fill- 
ing factor 5/2) with exact diagonalization utilizing both the spherical and torus geometries. The 
two distinct two-component systems we consider are the true bilayer and effective bilayers (wide- 
quantum- well). In each model (bilayer and wide-quantum- well) we completely take into account 
inter-layer tunneling and charge imbalancing terms. We find that in the half-filled lowest Landau 
level, the FQHE is described by the two-component Abelian Halperin 331 state which is remark- 
ably robust to charge imbalancing. In the half-filled second Landau, we find that the FQHE is 
likely described by the non- Abelian Moore-Read Pfaffian state which is also quite robust to charge 
imbalancing. Furthermore, we suggest the possibility of experimentally tuning from an Abelian to 
non-Abelian FQHE state in the second Landau level, and comment on recent experimental studies 
of FQHE in wide quantum well structures. 

PACS numbers: 73.43.-f, 71.10.Pm 



We theoretically investigate two-component spin- 
polarized fractional quantum Hall states in the bilayer 
half-filled first and second (first-excited) Landau levels 
considering the effects of both inter-layer tunneling and 
charge imbalance tunneling using exact diagonalization. 
The aspect that is different in this work compared to 
previous studies 1 " 5 is that we consider inter-layer tun- 
neling along with charge imbalance tunneling. The rea- 
son we consider this extra tunneling term (charge im- 
balance), along with the an inter-lay er t unneling term, 
is because recent experimental efforts^' have achieved 
bilayer FQHE systems where both the inter-layer and 
charge imbalance tunneling terms can be controlled by 
varying different system parameters (such as gate volt- 
ages) and a theoretical investigation of the physics of 
this system is both timely and important. Furthermore, 
the basic effects produced by a charge imbalance term 
in a theoretical exact diagonalization context is currently 
lacking. Before we delve into our results we provide an in- 
troduction, motivation, and historical perspective of this 
broad subject. 



I. FRACTIONAL QUANTUM HALL EFFECT IN 
THE ONE- AND TWO-COMPONENT VARIETY 

The fractional quantum Hall effect (FQHE) has 
proved to be the paradigm for emergent quantum 
physics for the nearly 30 years of its existence. It 
occurs when electrons are confined to a (quasi-)two- 
dimensional plane (such as in semiconductor structures, 



e.g., GaAs/AlGaAs, with electron densities on the or- 
der of 10 10 -10 11 cm -2 ) and a strong perpendicular mag- 
netic field is applied (usually on the order of tens of 
Tesla-sometimes up to ~ 40 T) . Phenomenologically, the 
FQHE manifests as a quantized plateau in the Hall re- 
sistance R xy (quantized to parts per billion) and a con- 
comitant vanishing (or deep minimum that displays acti- 
vated behavior) of the longitudinal resistance R xx . The 
FQHE is said to occur at rational fractional filling fac- 
tor v when the quantized value of the Hall resistance 
is R xy = h/{ye 2 ), where v = p/<fto (here p is the elec- 
tron density and <po = hc/eB is the magnetic flux quan- 
tum and B is the magnetic field strength, hence, v is the 
number of electrons per magnetic flux quanta) . The per- 
pendicular magnetic field quantizes the two-dimensional 
kinetic energy into Landau levels separated in energy by 
hu c = heB/mc and when the filling factor v is made to 
be fractional (like it is for the FQHE), by either adjust- 
ing the electron density and/or magnetic field strength, 
the kinetic energy is a constant and completely flat bands 
obtain. In the limit that huj c — > 00 (or the extreme quan- 
tum limit) the electron-electron Coulomb interaction is 
the dominant term in the Hamiltonian for electrons frac- 
tionally filling a Landau level (LL). 

The one-component FQHE in the lowest orbital elec- 
tronic Landau level is the most often discussed since 
it has been observed in the form of approximately 
80 odd-denominator FQHE states and is well under- 
stoooP^. Essentially, the FQHE occurs due to the non- 
perturbative and electron-electron interaction driven for- 
mation of an emergent topologically ordered^ incom- 
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pressible quantum fluid at certain filling factors v with 
non-zero energy gaps. This (bulk) energy gap, along with 
some sample disorder, explains the FQHEPh 

Recently, the FQHE in half-filled Landau levels has 
reinvigorated the community due to its possible con- 
nection to topological quantum computing, non-Abelian 
quasiparticles, and the requisite cutting edge material 
science advances that have produced much of this new 
physics and continue to coax nature into revealing her 
secrets. In particular, the FQHE at filling factor 5/2p3 
which corresponds to filling factor 1/2 in the second 
orbital electronic Landau level, has arguably produced 
most of the excitement. 

The denominator of v for the "standard" fractional 
quantum Hall states is always odd and reflects the 
fermionic nature of the quasiparticles. At z/=l/(even), 
naive, zeroth-order, theory asserts that no FQHE would 
be expected. This is because, at f=l/(even), the 
weakly inter acting quasiparticles of the FQHE (Compos- 
ite Fermion d 12 ! 13 *) experience a zero effective magnetic 
field and form a gapless exotic Fermi sea^HUD i n fact, 
no FQHE is experimentally observed in one-component 
systems at v = 1/2 instead showing signatures of the 
exotic Fermi seaPSEL] ( see Fig. [l]). 

The FQHE at v = 2 + 1/2 = 5/2 (the 2 comes from 
completely filling the spin-up and spin-down Landau lev- 
els) is particularly interesting and the only known and 
well established violator of the "odd-denominator" rule 
for the one-component FQHE. A one-component FQH 
state with an even-denominator suggests some kind of 
interaction driven pairing among the weakly interacting 
emergent (quasi-)fermions and does not find a descrip- 
tion within the standard FQHE theory. All FQHE states 
observed in the second Landau level (SLL), such as 5/2, 
7/3=2+1/3, 8/3=2+2/3, 12/5=2+2/5, etc. (there are 
about 8 FQHE states observed in the SLL in all) are 
fragile when compared to the FQHE states in the lowest 
Landau level. The 5/2 FQHE is one of the, if not the, 
strongest FQHE states in the SLL and yet has a measured 
activation gap only about 0.5 K or less. This is despite 
being measured in the world's cleanest two-dimensional 
electron gas samples with mobilities over 30 x 10 6 cm 2 /Vs 
at temperatures of less than 100 mK. Hence, there is 
a strong experimental drive to produce cleaner samples 
and apparatus capable of doing measurements at exceed- 
ingly low temperatures. Most early experimental obser- 
vations of the FQHE at 5 /2 indicated the state to be 
one-component in nature22H22] However, care needs to be 
taken when making broad abs olute st atements and there 
is some experimental evidence 1 22 * 27 * 30 ! that puts the one- 
component interpretation of the 5 /2 FQHE into question, 
but, the interpretation of these results is notoriously dif- 
ficultP. 

Recently, the FQHE was observed^ at v = 5/2 in a 
low magnetic field where the ratio of the cyclotron energy 
huj c to the Coulomb energy e 2 /el (where e is the electron 
charge, e is the dielectric constant of the host semicon- 
ductor and I = ch/eB is the characteristic length scale 
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FIG. 1. (Color online) (a) A Composite Fermion is an elec- 
tron bound (or attached) to an even number of quantum me- 
chanical vortices of the many-body wavefunction, see Ref. 1121 
and 1 131 (b) A compressible Fermi sea of Composite Fermions 
forms for a one-component system at filling factor v = 1/2 in 
the lowest Landau level and produces no FQHE since a (Com- 
posite Fermion) Fermi sea has no energy gap. (c) In the half- 
filled second Landau level at v = 5/2 the electron-electron 
interaction is modified compared to the electron-electron in- 
teraction in the lowest Landau level (see text) causing the 
weakly interacting Composite Fermions to pair into a spin- 
polarized p-wave BCS state described by the Moore-Read 
Pfaffian wavefunction. This state, due to the quasiparticle 
pairing, has an energy gap and thus exhibits the FQHE. 



called the magnetic length) is near unity or less. In that 
situation, it is not clear how the system would be com- 
pletely spin-polarized and/or would not experience sig- 
nificant Landau level mixing. Recent experiments22H251 
have begun to seriously tackle these two issues. Not sur- 
prisingly, theoretical groups have enthusiastically taken 
up the (considerable) challenge of understanding the role 
of Landau level mixmj^ jjjgl T hat being said, the FQHE 
at 5/2 has also been seer i 40 * 41 * at magnetic field strengths 
of more than 10 T. A FQH state at fields that high is 
likely to be spin-polarized and, hence, one-component. 
Addin g to the confusion is the fact that all theoretical 
wor ij42|43| £ Q date have established that, within a wide 
range of parameter space, the electrons are completely 
spin-polarized at v = 5/2. It is extremely difficult the- 
oretically (computationally) to consider fully spin un- 
polarized FQHE and until there is definitive experimen- 
tal evidence indicating the real system to be spin un- 
polarized we will theoretically consider the electron spin 
to be absolutely polarized throughout this work. 

Arguably, the reason the FQHE at 5/2 is so fasci nating 
is its connection to topological quantum computing 44 45 . 
An intriguing ansatz called the Moore- Read PfaffiarP^ is 
thought to describe the FQHE at 5/2 and this ansatz has 
non-Abelian quasiparticle and quasihole excitations. It is 
proposed^ that the world-lines of these non-Abelian ex- 
citations can be braided around each other, thus changing 
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the ground state in the degenerate manifold of ground 
states, and certain quantum computing gates can be 
achieved. This degenerate manifold of states is separated 
from the continuum by an energy gap that is topological 
in nature, since it is a FQHE state, thus, any sort of local 
disturbance to the system, like those caused from typical 
noise encountered in an experiment, will not be able to 
destroy the state due to its topological origin. Of course, 
if the disturbances are of an energy that is larger than the 
protective gap then the above does not hold. Using the 
non-Abelian quasiparticles and quasiholes of the Moore- 
Read (MR) Pfaffian (Pf) description of the 5/2 state of 
the FQHE to achieve quantum computing gates that are 
topologically protected (called fault-tolerant t opo logical 
quantum computation) is a major research goaPSl. (Note 
that the non-Abelian quasiparticles that might exist for 
the FQHE at 5/2 are so-called Ising anyons and are un- 
able to be used for universal quantum computation-see 
Ref. US] for more details on this point and how Ising 
anyons can be augmented to achieve a version of uni- 
versal fault-tolerant topological quantum computation.) 

Intuitively, the existence of the Moore-Read Pfaffian 
can be understood (see Fig. [I]). Any time a Landau level 
is half- filled (not just at v = 1/2 as mentioned above but 
also at v = 5/2 = 2 + 1/2) naive zeroth-order theory 
tells us that a (Composite Fermion) Fermi seaP^HUl can 
form. (Note that in this discussion we assume that the 
two completely filled Landau levels are inert which is the 
same as assuming that Tiuj c — > oo.) As is well known, any 
system of weakly interacting fermionic quasiparticles is 
unstable to pairing via the BC^m pairing mechanism 
and, if the quasiparticles are spin-polarized, the simplest 
pairing is chiral p-wave symmetry (p x +ipy) (see the work 
by Read and Green in Ref. 50 for a discussion of the 
BCS mean field description of the FQHE in a half-filled 
Landau level.) The resulting paired state has an energy 
gap and the FQHE will occur as long as the system is 
clean enough and cold enough such and the Fermi energy 
lies within the FQHE gap. The Moore-Read Pfaffian 
wavefunction encapsulates this physics. 

We b riefly note that recently it has been pointed 
out! 51 1 52 1 that a competing state for the FQHE at v = 5/2 
is the so-called anti-Pfaffian which is the particle-hole 
conjugate of the Moore-Read Pfaffian. The anti-Pfaffian 
is topologically distinct from the Pfaffian due to the fact 
that the two have different edge state behavior which, in 
principle, makes it possible to tell the two apart experi- 
mentally. We emphasize that both the Moore-Read Pfaf- 
fian and anti-Pfaffian both support non-Abelian quasi- 
hole and quasiparticle excitations so both could perform 
as platforms for fault-tolerate topological quantum com- 
putation. Whichever one happens to be responsible (if 
either) for the FQHE at v = 5/2 depends on how, and if, 
the particle-hole symmetry is broker^ as well as experi- 
mental and material details such as disorder and Landau 
level mixing^2ll. 

(All results in this work concerning the Pfaffian apply 
equivalently to the anti-Pfaffian within the constraints of 



our approximation scheme since the model we use does 
not distinguish between Pfaffian and anti-Pfaffian. The 
reason for this is that our Hamiltonian is particle-hole 
symmetric, since we do not consider any particle-hole 
breaking terms such as those that might arise due to 
Landau level mixing. Since the Pfaffian and anti-Pfaffian 
are particle- hole conjugates of one another, their (bulk) 
physics, the subject of our current work, is identical.) 

A natural question is why the FQHE is observed at 
v = 5/2 = (2 + 1/2) and not, so far, at v = 1/2 for 
one-component systems. Theoretically, this is due to the 
physics of the lowest Landau level being different com- 
pared to the physics of the second Laudau leveS^HI3 i. e ., 
details matter. For the FQHE at filling factor 5/2, the 
inert electrons in the lowest spin-up and spin-down Lan- 
dau levels partially screen the electron-electron Coulomb 
interaction between the electrons half-filling the second 
Landau level-in fact, the electrons in the lower Landau 
level over screen the Coulomb interaction 56 which leads 
to a slight attraction among the quasiparticles causing 
them to form a BCS state-the Moore- Read Pfaffian state. 
Furthermore, the experimental details of a real quan- 
tum well-the single-electron wavefunction in the direc- 
tion perpendicular to the two-dimensional surface has a 
finite extent (the so-called finite thickness of the quasi- 
two-dimensional electron system where typical experi- 
mental systems have widths ranging from approximately 
~ 20 nm for the thinnest samples to ~ 60 nm for the 
wide-quantum-well samples) which produces subtle ef- 
fects that make the effective interaction felt by the elec- 
trons in the half-filled second Landau level even more 
amenable to forming a non-Abelian Moore-Read Pfaf- 
fian FQHE at v — 5/2. Recently, it has been theoret- 
ically shown that (within certain approximate models) 
the effects of Landau level mixing^tlH might produce 
non-trivial and subtle effects that drive the system to 
the PfaffiarPl or the anti-PfafnaiP 1 state. 

We reiterate that experimentally the FQHE at v = 1/2 
has n ot b een observed in one- component systems and 
theorjBEsl suggests that it will not occur for a pure 
Coulomb interaction. However, one cannot rule out a 
Moore-Read Pfaffian FQHE at v = 1/2 if the system pa- 
rameters are tuned in the perfect waySSH6i] or) p er haps, at 
extremely low temperatures not currently experimentally 
accessible. 

However, about the same time the experimental obser- 
vation of the 5/2 FQHE was realized, the FQHE was also 
observed by Suen et alW^^ and Eisenstein et alW^ in the 
lowest Landau level at filling factor v = 1/2 in systems 
that were later determined by He et aZP"^ to be spin- 
polarized two-component systems. It turns out that the 
FQHE at v = 1/2 in two-component systems is described 
by the (Abelian) Halperin 331 two-component wavefunc- 
tiorP3. The system of EisensteirP^ was an actual bilayer 
(two quasi-two-dimensional systems separated from one 
another by a tunneling barrier-typical parameters are 
quasi-two-dimensional electron systems of width ~ 20 
nm separated from one another by a tunneling barrier of 
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thickness ~ 3 — 10 nm with electron densities on the or- 
der of 10 11 cm -2 ), thus, the two-components are the two 
layers. The FQHE at v — 1/2 in bilaycrs occurs by a two- 
part process where LaughlirP states, at filling factor 1/3, 
are formed in each layer and the fermions between lay- 
ers form pairs. In the experiments by Suen et aZP^HsU the 
electrons in the wide-quantum- well (a width of order ~70 
nm and electron densities of ~ 10 11 cm -2 ) minimize their 
single-particle energy by reorganizing into effectively two 
"layers" and, again, the FQHE at v = 1 /2 occurs as de- 
scribed above. Recent observations^^ of v = 1/2 and 
1 /4 in wide-quantum- wells have rekindled interest in this 
rich system. 

The discovery and identification of the FQHE at v = 
1/2 was an exciting advance in the FQHE because it 
opened up the possibility of a new fractional quantum 
Ha ll state o utside those given in the "standard" the- 
Or ^9 | io | i2 ji3] However, the Halperin 331 state is closely 
related to the standard theory being that it is essentially 
a pairing of Laughlin states. Scarola and Jair]2£l further 
generalized the Halperin 331 state to pairings of FQHE 
states belonging to the Jain sequence producing states 
for bilayer systems at filling factors other than v = 1/2- 
they also produced partially pseudo-spin polarized states 
at v = 1/2 that are distinct from the Halperin 331 state. 

Along with the recent observations^^ of v = 1/2 and 
1 /4 in wide-quantum- wells are observations^^ of bilayer 
FQHE at filling factor 1/2 in the lowest Landau level 
with asymmetric charge distributions-so-called "tilted 
samples". A tilted sample is created by an additional 
charge gate that produces a charge asymmetry by "push- 
ing" electrons from one side of the wide-quantum- well to 
the other. In the simplest approximation, this asymme- 
try produces a charge imbalance between one layer and 
the other while keeping the total filling factor v fixed. 
That is, the individual filing factors in each layer, and 
i^2, respectively, (keeping v\ + v 2 = v fixed) can be var- 
ied between Lvi,!^) = (^,0), to {v\,v%) = (v/2,v/2), to 

Of course, care should be taken when considering the 
effects of a charge gate that presumably causes charge 
asymmetry in a wide-quantum-well. Recently, Scarola et 
aZP^have used the local-density-approximation (LDA) to 
get a more realistic handle on exactly what tilting does to 
the charge distribution and relative single-electron energy 
levels in the quantum well. They went on to analyze the 
system with the use of variational wavefunctions finding 
that the recent experimental observation of the FQHE 
at v — 1/2 by Shabani et al. 7 at intermediate charge 
imbalancing is likely describ ed by a partially psuedo-spin 
polarized bilayer stat d 68 | 69 J . 

The purpose of the present work is to understand the 
general effects of charge imbalancing in the most minimal 
model one can consider to describe the bilayer FQHE 
in an exact context-that is, using exact diagonalization. 
Our work compliments recent workP^ since we are solving 
the Hamiltonian exactly instead of using a combination 
of LDA and variational wavefunctions. As such, we are 



restricted to only certain ansatz, namely the Moore-Read 
Pfaffian and the Halperin 331 state. We emphasize that 
our solution is general and exact. 

In Sec. [H] we present our theoretical model in the form 
of a Hamiltonian that can describe either a true bilayer 
and effective (wide-quantum-wcll) bilayer where we con- 
sider both inter-layer and charge imbalance tunneling 
terms. Four (two plus two) variational states are dis- 
cussed in Sec. |III| that are thought to describe our bi- 
layer and effective bilayer Hamiltonian. Sees. |IV| and [V] 
present our results for the bilayer and wide-quantum- well 
model in both the lowest Landau level (in both the spher- 
ical and torus geometry) and second Landau level (in the 
spherical geometry). Finally, in Sec. VI we present our 
conclusions. 



II. THEORETICAL MODELS: TRUE BILAYER 
AND EFFECTIVE BILAYER 

There are two experimental systems in which two- 
component FQHE systems can be produced. One of them 
is to manufacture a true bilayer system which consists of 
two parallel (quasi-)two-dimensional electron systems of 
width w separated from one another by a tunneling bar- 
rier of thickness d. This system is used extensively by 
the Eisenstein group at CalTecfPl and is what one gen- 
erally thinks about when contemplating a bilayer system. 
The height of the barrier can be adjusted such that elec- 
trons are either localized in separate layers or delocal- 
ized between the two layers, i.e., large energetic barrier 
or small energetic barrier, respectively. In other words, 
the symmetric-anti-symmetric energy gap, As^g, can be 
controlled independently of the individual well width w 
and the layer separation d. In this system, an electron 
can be in the right (R) or left (L) state, written as \R) or 
\L), respectively. Fig. [2] shows a schematic diagram of a 
bilayer system consisting of two parallel two-dimensional 
electron systems where the electrons in each layer are 
confined to two-dimensions by a quantum well of width w 
(both layers have the same width) and the quantum-well- 
center to quantum- well-center separation is d (d > uu/2). 
A typical density profile of electrons in the right or left 
layer ((R\R) or (L\L)) is also shown in Fig. [2^. The 
ground state of the single-electron bilayer Hamiltonian is 
the symmetric state S given by 



|S) = 



\R) + \L) 
V2 



(1) 



and the next higher energy state is the anti-symmetric 
state given by 



\AS) = 



\R) ~ \L) 
V2 



(2) 



(assuming no charge imbalance term in the Hamiltonian). 
The two states (S and AS) are separated by an energy 
difference Ag^s and their respective (typical) density 
profiles ((SIS) and (AS\AS)) are shown in Fig. [2b. 
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FIG. 2. (Color online) A bilayer system consisting of two 
(quasi-)two-dimensional quantum wells of width w separated 
from one another by d (d > w/2). (a) Typical density profile 
of an electron localized in either the right (R) or left (L) layer, 
respectively, (b) The ground state single particle wavefunc- 
tion (the symmetric (S)) state and the next higher energy 
state (the anti-symmetric (AS)) state separated in energy 
from one another by the so-called symmetric-anti-symmetric 
energy gap A S as- 



The other way to experimentally create a two- 
component FQHE system is to make a single (quasi-)two- 
dimensional electron gas of large width, a so-called wide- 
quantum- wellTWQW), as done by the Shayegan group 
at Princeton 6 7 62 64 . In this case, the electrostatic inter- 



action between the electrons (at the Hartree-Fock level) 
compels the system to behave as an effective bilayer. 
The electron density for the ground state is maximum 
(and symmetric) near the edges of the wide-quantum- 
well and depleted in the middle. The energy level di- 
agram and typical density profiles are similar to those 
shown in Fig. |3^. Note that in this system, As^g is a 
strong function of the width W (and electron density) of 
the wide-quantum-well. This contrasts the bilayer sys- 
tem described above where, in principle, As^s, the in- 
dividual well width w, and the layer separation d can all 
be independently adjusted. 

A simplified model for the WQW is due to Papic et 
aZP where the ground state of a wide-quantum-well of 
width W is taken to be a symmetric state (for example, 
(z\S) ~ sui(itz/W)) and the next excited state an anti- 
symmetric state (for example, (z\AS) ~ sm(2wz/W))- 
both the S and AS states are written on the inter- 
val z 6 [0, W\. (The coordinates of the (quasi-)two- 
dimensional plane will always be denoted by x and y 
and the coordinates of the wide-quantum-well width, bi- 
layer layer separation, or individual quantum-well thick- 
ness will be denoted by z.) The density profiles for the 
WQW model are shown in Fig. [3]d. In this model, one 
can transform to the layer-basis (right and left layers) by 
taking \L) = (\S) + \AS))/V2 and \R) = (\S)-\AS))/y/2, 
cf. Fig. [J 

The Hamiltonian for our two-component model sys- 
tems can be written (in the S'AS'-basis) as 



H 



^ (™iCTi,m2CT2|^|m 3 CT3,m 4 CT4)c^ iCTi c^ 2(T2 c m 3 (T3 c m4Cr4 (S mi+m ^ m3+m4 

{rm,<ri=S,AS} 

^SAS sr^f t t \ , ST^i t , t \ 

~Y~ l^( C ™S Cm S ~ c' mAS C mAS ) + — }_^(c ] mS C m AS + c' mAS C mS ) . 



(3) 



The second to last term (the term with the prefactor Asas/%) controls inter-layer tunneling and, in the S'AS'-basis 
is represented by the pseudo-spin operator S z , while the last term (prefactor A.p/2) controls charge imbalance and is 
represented by the pseudo-spin operator S x . In our convention, positive Ap drives the electrons into the right (R) 
layer. Notice that (mi<7i, m2(J2\V\m^a^, m^a^ is different in the bilayer case than it is in the WQW case, in general. 
Namely, in the bilayer case, the Coulomb matrix element is found by using the transformations given in Eqs.[T]and[2] 
and the fact that the potential energy between two electrons a distance r apart in a given layer (the right layer, say) of 
width w is given by the Zhang-Das Sarma 7 - potential V5 ntra (r) = 1/\A" 2 + w 2 and the potential between two electrons 
in two different layers separated by a distance d is given by ^ n t e r(f) = 1/vV 2 + d 2 . In the WQW system the Coulomb 
matrix element is calculated by considering the wavefunctions for the S and AS states ((z\S) and (z\AS)) and doing 
the required integrals. 

In the (left/right) layer-basis we can write the bilayer Hamiltonian very simply as 



H = Y^i V intU\r l - Ijl) + VtatraG?* ~ ^|) + Knterd^ - Tj\)] - A SAS S X + A P S Z 



(4) 



where operators O with a tilde are written in the layer-basis and coordinates r and f belong to electrons in different 
layers. Switching from the layer- to the S'AS'-basis is a pseudo-spin rotation. Note that for d = the Hamiltonian is 
symmetric between As^s and Ap. Increasing d destroys this symmetry and it becomes harder to drive the system 
to be one-component in the layer-sense. 

I 



In all of our results we use the natural FQHE units: lengths are given in units of the magnetic length I and 
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FIG. 3. (Color online) (a) A two-component system created by a wide-quantum- well (WQW). The electrostatic potential 
generated by the electrons self-consistently creates an effective bilayer system by pushing the electrons towards the walls of the 
WQW. Note that (a) is not identical to Fig. [2ja even though they are very similar-specifically, As as is a strong function of W 
and the electron density in WQW systems as opposed to bilayers. (b) The model of Papic et aZ™ which can be transformed 
into the layer-basis shown in (c). 



energies are given in units of the Coulomb energy e 2 /el. 

(Note that the system we are considering is a rather 
general two-component Hamiltonian and when the layer 
separation is zero any results we find are applicable to 
SU(2) symmetric two-component systems where the con- 
stituents interact by a Coulomb-like interaction-exactly 
Coulombic for the lowest Landau level and a slightly 
modified Coulombic interaction for the second Landau 
level. For example, when d = the Hamiltonian (Eq. [3| 
describes spin-full electrons interacting via a Coulomb 
interaction with both S z and S x terms.) 

III. VARIATIONAL WAVEFUNCTIONS: 
MOORE-READ PFAFFIAN AND HALPERIN 331 

We consider four variational wavefunctions to describe 
our FQH system for half-filled lowest and second Lan- 
dau levels {y — 1/2 and v = 5/2, respectively) even 
though it will appear at first that we are only consider- 
ing two. Namely, we consider the one-component Moore- 
Read PfaffiarpEl (Pf) and two-component Halperin 331 
(331) wavefunctions^. (Note that worlPS^ has been 
done on the FQHE at total filling factor of unity, i.e., 
i/ = 1/2 + 1/2=1, however, that is not the case that we 
study in this work.) 

The Moore-Read Pfaffian state is written as 

f 1 1 N 
*Pf - Pf = \ i[(z t ,a - z^ a f (5) 

where z^ a = Xi. a — iyi, a is t he position of the ith electron 
in complex coordinates with a labeling its state, i.e., a = 
S 1 , AS 1 , R or L. The origin of this wavefunction can be 
understood at an intuitive level. If one writes down an 
ansatz wavefunction to describe a BCS paired state of 
Composite Fermions in a half-filled Landau level then 
one arrives at a wavefunction of the Pfaffian forrrF^, i.e., 

Pf{s(r«)}n(*-*i) a = 

i<j 



-^{3(^12)3(^34) • ■ ■ g(rN-i,N)} Y[( z i ~~ z ^ 2 ( 6 ) 

i<j 

where g(\vi — rj\) is the real space pairing amplitude 
for a pair of electrons located at position and r^, 
rij<j( z j ~ z j) 2 binds two quantum mechanical vortices of 
the many-body wavefunction to each electron (the Com- 
posite Fermion transformation and fixes the filling factor 
in a Landau level to 1/2). A is the anti-symmetrization 
operator. 

Despite the intuitive picture that goes along with the 
Moore-Read Pfaffian wavefunction it was originally de- 
rived using conformal field theory and, as such, the Pf 
has a low-energy effective conformal field theory descrip- 
tion, cf. Refs. Bol, 1751 and I5U1 It is within this conformal 
field theory that the understanding of the non-Abelian 
nature of the quasiparticle excitations was first illumi- 
nated. However, whether or not the Pf has anything to 
do with the physics of the FQHE at 5/2 has depended 
crucially on numerical calculations and compari sons with 
the results of exact diagonalizatio n 42 * 43 l 58 " 61 l 80 l As men- 
tioned above, the details matter greatly in determining 
the p hysics of the FQHE at 5/2, such as, finite thickness 
effect s 59 ! 60 !, the competition between FQHE states and 
non-FQHE striped phases ^ an d the compressible (Com- 
posite Fermion) Fermi se a 1 16 ! 17 ! (as mentioned, the latter 
occurs in the lowest Landau level at v = 1/2 while the 
former occurs in the second Landau level at v = 5/2), 
the effects of Landau level mixinj2SH2U ; e t c Note that 
most of the numerical studies mentioned above do not 
treat the v = 5/2 problem in full complexity, i.e. with 
the lowest Landau level filled with two spin species and 
one Landau level half filled. Instead, because of computa- 
tional simplicity, one considers a half filled lowest Landau 
level (usually without spin) with an effective interaction 
that contains the form factors of the first excited Landau 
level-the appropriate Haldane pseudopotentials 81 . This 
is what we mean by LL v = 5/2 in an exact diagonalization 
context. 

Thus, it should be clear that any effective conformal 
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field theory (or effective BCS mean field HamiltoniarP^) 
written down to describe the physics at 5/2 is only as 
good as its physical predictions and agreement with ex- 
periments and its agreement with numerical calculations. 
This is because the entire reason for the existence of the 
FQHE at all filling factors is due to non-perturbative 
physics arising from the strongly interacting electrons in- 
teracting with a Coulomb interaction modified by the de- 
tails of which Landau level is fractionally filled by elec- 
trons, the thickness the quantum well, the amount of 
Landau level mixing that is taking place, etc. Hence, 
any mean field theory that throws away all these details 
at the first step cannot explain, for example, why the 
FQHE occurs at v = 5/2 and not v — 1/2 (or for that 
matter, at 9/2 or 13/2...) without the aid of both numer- 
ical calculations and, most importantly, experiments. In 
the end, whether a particular candidate wavefunction, 
e.g. Moore-Read Pfaffian or Laughlin or Halperin 331 or 
Composite Fermion states, etc., describes a real FQHE 
at some filling factor is an energetic question delicately 
depending on competing states where all the details of 
the microscopic Hamiltonian may eventually matter, and 
therefore, extreme caution is necessary in identifying ex- 
perimental FQHE states with candidate incompressible 
wavef unctions, particularly in higher LLs and/or multi- 
component systems where various competing states are, 
in general, more viable. 

The Halperin (two-component) 331 wavefunction is 
written as 



(a) 



(b) 



leftlayeX^ ^^right layer 
•Ppf (layer-basis) 




*Ppf (SAS-basis) 



(C) 



L 

4t 



(d) 




^¥iM (layer-basis) 



"P 3 3i (SAS-basis) 



*331 = 

N/2 N/2 N 

Y[(zi,i - z jt i) 3 \[ i ~ -' If - 2i,a) , (7) 

i<j i<j i,j 

where the subscript 1 (or 2) on z^x labels the quantum 
number of the electron. This quantum number could de- 
scribe spin (up or down), layer (R or L), subband (S 
or AS). Note that in our work, we consider the bilayer 
situation so the subscript will either denote the layer or 
the symmetric or anti-symmetric state depending on in 
which basis we choose to work. Intuitively, the origin of 
this wavefunction is that, in a bilayer system in a half- 
filled Landau level, the electrons form v = 1/3 Laughlin 
states^ dli<j( z i — z j) 3 ) am ong each electron component 
and then pair among different components (simple Jas- 
trow factor H^z^i - z j>2 )). 

(Note that we always drop the Gaussian factors (~ 
exp(— | -^i , a, | 2 ) ) that are always present when writing 
wavefunctions describing electrons entirely in (or pro- 
jected into) the lowest Landau level-in fact, the Gaus- 
sian is often considered to be part of the measure so one 
is technically not "dropping" anything at all.) 

It is important to realize that the Moore-Read Pfaf- 
fian is a one-component state that can be one-component 
in both the S'AS'-basis sense or in the layer-basis sense, 
see Fig. [4^ and |4Jd. That is, the Pf wavefunction either 
describes pairing among one-component electrons in the 



FIG. 4. (Color online) The Moore-Read Pfaffian wavefunc- 
tion written in the (a) layer- and (b) SAS-basis where the 
p-wave pairing in the former is between Composite Fermions 
in the right (or left) layer and the p-wave pairing in the SAS- 
basis is between electrons in the symmetric (S) state. The 
Halperin 331 state is shown in (c) the layer-basis and (d) the 
SAS-basis where in the layer-basis the 331 state pairs Com- 
posite Fermions between layers and in the SVlS-basis it pairs 
Composite Fermions in the symmetric (S) and anti-symmetric 
(AS) states. 



right (R) layer or left (L) layer (the layer-basis) or pair- 
ing among one-component electrons in the symmetric (S) 
state (the SAS'-basis). The former is the Pfaffian state 
that has its origin in the single-layer FQHE, while the 
latter is the Pfaffian understood as the even (symmetric) 
channel of the BCS description by Read and GreerPSl. 
(We choose our charge imbalancing term to drive the 
electrons into the right layer, however, we could have 
equivalently chosen the charge imbalance term to drive 
the electrons into the left layer.) Thus, there are actually 
two Moore-Read Pfaffian states to consider when both 
tunneling terms (Asas and Ap) are non-zero (or three if 
we allow the sign of Ap to change). When the system is 
SU(2) symmetric, i.e., when the layer separation is zero, 
the two different states are simply related via a rotation 
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in pseudo-spin space. 

The Halperin 331 state can also be two-component in 
two ways: it can describe pairing of electrons in the R 
and L layers (the layer-basis) or it can describe pairing 
of electrons in the S and AS states (the SAS'-basis), see 
Fig. [4p and Fig. [4]d. Hence, there are two Halperin 331 
states to consider when both tunneling terms are finite. 
Note, however, that the Halperin 331 state, as usually 
understood, is defined in the layer-basis. Again, in the 
SU(2) situation, these two 331 states are related via a 
pseudo-spin rotation. 

Thus, instead of dealing with two variational wave- 
functions (Pf and 331) we are dealing with four-two Pfs 
and two 331s. Even though the two pairs of wavefunc- 
tions are related via a pseudo-spin rotation, when the 
SU(2) symmetry is broken this is no longer the case and 
the determination of which state describes the physics is 
non-trivial. 

In our calculations we utilize the spherical geome- 
try^ where the electrons are confined to the surface of 
a sphere of radius J N^/2 and a radial magnetic field 
(perpendicular to the surface) is produced by a magnetic 
monopole of strength N^/2 placed at the center. (We 
briefly consider the torus geometry below.) The total 
magnetic flux piercing the surface is and is either an 
integer or half-integer according to DiracP^l The filling 
factor (in the partially filled Landau level) is given by 
liniAr_ i . 00 N/Nff, for N electrons. In our case, we are con- 
sidering the Moore- Read Pfaffian and Halperin 331 states 
which describe a half-filled Landau level, thus, the rela- 
tionship between N and for a finite spherical system 
is Ntf, — 2N — 3. Due to computational constraints we 
consider N = 8 and N$ = 13 (the N = 6 particle sys- 
tem is aliased with a 2/3 filled Landau level and, thus, 
produces ambiguous results and the N = 10 electron 
system is too big to consider while adequately exploring 
the large parameter space inherent when tackling a two- 
component Hamiltonian with both inter-layer and charge 
imbalancing tunneling). In the spherical geometry, pos- 
sible FQHE states are uniform states with total angular 
momentum L = 0. 

In this work we consider the overlap between the ex- 
act ground state of the model Hamiltonian and the four 
variational states: (1) Pf in the S'AS'-basis, (2) Pf in the 
layer-basis, (3) 331 in the SAS-basis and (4) 331 in the 
layer-basis. 



IV. RESULTS: LOWEST LANDAU LEVEL 

We first present our results for the lowest Landau level 
where we consider (i) the wavefunction overlap between 
the four variational states (Pf in the layer- and SAS- 
basis and the 331 in the layer- and STlS-basis) and the 
exact ground state of the Hamiltonian-note that since 
all the variational states describe the FQHE they have 
L = and so if the exact ground state of the Hamilto- 
nian does not also have total angular momentum L = 



the overlaps will be trivially zero since they are then of 
different symmetry, (ii) the pseudo-spin expectation val- 
ues of the exact ground state, which is more descriptively 
referred to as the expectation value of (Ns — Nas)/2 or 
(N R -N L )/2 where N s , N AS , N R and N L are the expec- 
tation values of the number of electrons in the symmetric 
state, the anti-symmetric state, the right layer and the 
left layer, respectively. Alternatively, this is simply (S z ) 
and (S x ) in the ^A^-basis, but we prefer to use the more 
physical (N s - N AS )/2 and (N R -N L )/2 since this does 
not depend on the basis in which we choose to write the 
Hamiltonian. (iii) the FQHE energy gap which is defined 
here as the energy difference between the L = ground 
state and the first excited state (this is often times re- 
ferred to as the "neutral gap" but we will call it the 
FQHE gap throughout this work) and, lastly, (iv) we in- 
vestigate the energy spectra in the torus geometry. We 
show all of these quantities from a variety of vantage 
points to give the clearest picture of the physics since 
things can get complicated rather quickly with so many 
parameters in the Hamiltonian, namely, layer separation 
d or WQW width W, SAS energy gap Asas and the 
charge imbalancing gap Ap. 



A. Bilayer 

We first present our results for the bilayer Hamilto- 
nian. Fig. [5^, shows the overlap between the Moore-Read 
Pfaffian and the exact ground state in the layer-basis (left 
column) and S'AS'-basis (right column) as a function of 
Asas and Ap for values of layer separation d ranging 
from d = 0.05 to 6. For the layer-basis (left column), it 
is clear that increasing Ag^g does not drive the system 
into the Pf phase, i.e., it does not increase the overlap 
between the layer-basis Pf with the exact ground state. 
However, increasing the charge imbalance Ap does drive 
the system into the Pf phase. For increasing values of d it 
takes a larger value of Ap to produce a ground state with 
a high overlap with the layer-basis Pf. This behavior is 
easy to understand. Non-zero Ap makes the system one- 
component in the layer-sense and thus a layer-basis Pf 
wavefunction (which is a one-component wavefunction) 
has a large overlap. It is also understandable that finite 
layer separation d would make it harder to polarize the 
electrons in the layer-sense, i.e., the electrons prefer to re- 
main two-component so they can take advantage of the 
reduction in potential energy between electrons in neigh- 
boring layers (the interaction is 1/ \Jr 2 + d 2 as opposed 
to 1/r for electrons in the same layer). 

The overlap between the exact ground state and the 
Moore-Read Pfaffian in the S^S-basis shows essentially 
the opposite behavior. This time, non-zero As^s drives 
the system to be one-component. Notice that for the 
smallest value of layer separation d = 0.05 the system 
is nearly SU(2) symmetric and the two bases (layer and 
SAS) are (nearly) pseudo-spin rotations of each other. 
The figure is deceptive because the scale on the two axes 





FIG. 5. (Color online). Lowest Landau level: (a) Wavefunction overlap between the exact ground state of the bilayer Hamil- 
tonian and the Moore-Read Pfaffian written in the layer-basis (left column) and the .SMS-basis (right column) shown as a 
function of inter-layer tunneling Asas and charge imbalance Ap for different values of layer separation d and zero individual 
layer thickness w = 0. (b) Same as (a) but for the Halperin 331 wavefunction. (c) Pseudo-spin expectation value or, more 
physically, the expectation value of (Nr — Nl)/2 (left column) and (Ns — Nas)/2 (right column) as a function of Asas and 
Ap. (d) shows the FQHE (neutral) energy gap for the exact bilayer Hamiltonian. 



is different, but the two figures are essentially identical 
upon reflection across the AgAS = Ap line. This sym- 
metry is quite obviously destroyed upon increasing d. 

A naive look at Fig. [5^, would lead to a conclusion 
that we have found the Pfaffian state in some regions 
of Asas — Ap phase diagram. The reader may wonder 
how robust these conclusions are with increasing system 
size. Although we are unable to calculate the full phase 



diagram for A = 10 electrons, we can gain insight into 
the limiting cases when either Ap or Ag^g is very large. 
If we work in the layer-basis and Ap is very large, we 
recover the single-layerphysics; for our system of A = 8 
electrons, it turns outpQ that the ground state for the 
lowest Landau level Coulomb interaction is already in 
the universality class of Moore-Read Pfaffian. However, 
this no longer holds for larger systems such as A = 10 
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FIG. 6. (Color online) Lowest Landau level: Wavefunction 
overlap between the Moore-Read Pfaffian wavefunction in the 
layer-basis (top panel) and the S'AS-basis (lower panel) and 
the exact ground state for (a) d — 0, (b) d = 1, (c) d = 2 and 
(d) d = 5. 





FIG. 8. (Color online) Lowest Landau level: (Pseudo-spin) 
expectation value of the exact ground state of (Nr — N^)/2 
(top panel) and (Ns — Nas)/2 (lower panel) for (a) d = 0, 
(b) d = 1, (c) d = 2 and (d) d = 5. Note that these figures 
(especially for d = 0) are qualitatively similar to our cartoon 
schematic in Fig. [To] since the psuedo-spin expectation value 
essentially describes the one- or two-component nature of the 
ground state. Also, finite d breaks the SU(2) symmetry of 
reflection across the Asas = Ap line assumed in Fig. |10| 



(a) LLL, rf=0 (b) LLL, d=\ 




FIG. 9. (Color online) Lowest Landau level: FQHE energy 
gap for (a) d = 0, (b) d = 1, (c) d = 2 and (d) d = 5. 



FIG. 7. (Color online) Lowest Landau level: Wavefunction 
overlap between the Halperin 331 wavefunction in the layer- 
basis (top panel) and the SAS'-basis (lower panel) and the 
exact ground state for (a) d = 0, (b) d = 1, (c) d = 2 and (d) 
d = 5. 



and N = 12, leading to the conclusion that high over- 
laps with the Pfaffian in the layer-basis in Fig. [5^ are 
a finite size effect for v = 1/2. On the other hand, if 
we work in the S'AS'-basis and increase As^s, all the 
electrons will occupy the symmetric subband. This is 
again a one-component system, but with the interaction 
slightly softened with respect to pure CoulomtP 5 . It is 
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A-SAS " 



FIG. 10. (Color online) By varying either Asas or Ap the 
system can be two-component (in the SAS'-basis or layer- 
basis sense) when Asas and Ap are small and one-component 
(again in either the SVlS-basis or layer-basis sense) for large 
values of either Asas or Ap. The shading qualitatively indi- 
cates the one-component to two-component nature of the sys- 
tem with the darker shade indicating a two-component system 
and the lighter color a one-component system. 



knowrF^ that softening the Coulomb interaction via vari- 
ous finite- width corrections can lead to a phase transition 
between compressible and incompressible (non-Abelian) 
states, but this also leads to a significant reduction of the 
gap (as we show below). Therefore, the Pfaffian which 
is seen in our data for SAS basis is also not a typical 
incompressible state, but rather a "critical" state with 
some Pfaffian correlations and a small gap. In the re- 
mainder of this section, we will nonetheless, for the sake 
of brevity, refer to the one-component phase (fully polar- 
ized by either Asas or Ap) as described by the appro- 
priate basis "Pfaffian" state, cautioning the reader that 
such a designation most likely does not survive in the 
thermodynamic limit. 

Next we consider the overlap between the exact ground 
state and the Halperin 331 state in the layer-basis (left 
column) and 5^4<S-basis (right column) in Fig. [5}d. For 
the S'AS'-basis 331 state the overlap is very small for all 
values of AsaSi Ap or layer separation d. On the other 
hand, the overlap with the layer-basis 331 state is non- 
zero for moderate values of A 5 as and for 1 < d < 4 the 
overlap is large and, despite the fact that non-zero Ap 
eventually destroys the 331, it is remarkably robust to 
charge imbalancing. 

In Figs. [6] and [7] we show a clearer view of the overlaps. 
In particular, for d = when the system is SU(2) sym- 
metric it is obvious in Figs. [6] and [7] that the layer-basis 
and S^lS'-basis are merely pseudo-spin rotations of one 
another. Furthermore, it is clear that for d = there is 
no region in phase-space were the Halperin 331 wavefunc- 
tion (in either basis) is a good description and only upon 
increasing d away from zero do we see any sort of rea- 
sonably sized overlaps when increasing Asas (though it 
should be noted that for d = 1 and 2 the overlap with 331 



is sizable for small values of both Asas and Ap). For the 
half-filled Landau level at d = 0. th eoretically, the system 
is known to be spin-polarizecP^HSl, or one- component, so 
it is not surprising that a two-component wavefunction 
like the Halperin 331 state would not be a good descrip- 
tion and, therefore, have a small overlap with the exact 
ground state. Note also that the one-component Moore- 
Read Pfaffian is also not a good description for the d = 
case since, in the lowest Landau level, the ground state 
is most likely a non-FQHE (Composite Fermion) Fermi 
sea™]. 

Generally, on the example of our N — 8 system we 
can conclude that (i) when the system is largely two- 
component in nature then the overlap between the exact 
ground state and the Halperin 331 state is large, and 
(ii) when the system is largely one-component in nature 
the overlap between the exact ground state and the Pf 
is large. This is more clearly shown in Fig. [5]; where in 
the left and right columns we show the expectation value 
of (N R - N L )/2 and (N s - N AS )/2, respectively. When 
the system is largely one-component, either in the layer 
sense or SAS sense, the system concomitantly has a siz- 
able overlap with the one-component Pf state and when 
the system is largely two-component the overlap with 
the Halperin 331 state is large. Again, Fig. [8] shows the 
pseudo-spin expectation values clearer for a few specific 
examples of layer separation <i-note again how non-zero d 
destroys the symmetry between the layer- and S.AS-basis 
obtained when d = and the system is SU(2) symmet- 
ric. We emphasize that the first conclusion (i) is not 
affected by finite-size effects, whereas the conclusion (ii) 
likely holds only for a few special systems such as N = 8 
and it may be a finite-size artifact. 

In Fig. [TO] we show a schematic diagram that encap- 
sulates broad features of the bilayer model (and WQW 
model) for small layer separation d (small WQW width 
W). For large Asas and small Ap the system is one- 
component in the S'A^'-basis and for large Ap and small 
Asas the system is one-component in the layer-basis. 
When, both tunneling strengths approach zero, the sys- 
tem is two-component in the layer-basis, the 5j45-basis, 
or both. The diagram would be topologically similar for 
non-zero (and even large) values of layer separation d and 
WQW width W . The only difference is that the diagram 
would not be symmetric between Asas and Ap since the 
system is not SU(2) symmetric. 

Lastly, in Fig. [5ji we show the FQHE energy gap. For 
d ^$ 1 the FQHE energy gap is relatively constant as 
Asas and Ap vary. When d is increased past approxi- 
mately d ~ 2 the FQHE gap shows interesting features 
as a function of the strength of the two tunneling terms. 
Somewhat surprisingly, the FQHE gap obtains a maxi- 
mum for non-zero Asas and finite d (this is related to 
similar results found recently in Ref. [5] where the FQHE 
gap is maximum on the "ridge" as a function of Asas 
and separation d for zero charge imbalance Ap). This 
"ridge" basically separates the regions in the quantum 
phase diagram where the system is either in the Pf or 
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FIG. 11. (Color online) Lowest Landau level: Same as Fig. [5] except that all plots are shown as a function of inter-layer 
tunneling Asas and layer separation d for a number of different charge imbalances Ap and zero individual layer thickness 
w = 0. 



331 phase. Fig. [9] is a more detailed depiction of the 
FQHE energy gap. 

Before moving on to results for the second Landau level 
we show essentially the same results as Fig. [5^,-d but as a 
function of Ag^s and layer separation d for several values 
of Ap in Fig. [TTfc-d. This presentation is to more easily 
compare with previous bilayer worka^ where the data 
was presented this way. When considering the overlap 
between the exact ground state and the Pf we see that in- 
creasing imbalance (increasing Ap) eventually produces 
a high overlap for small values of d in the layer-basis. 
In the S'AS'-basis, charge imbalance pushes the maxi- 



mum overlap to larger values of Ag^s and slightly larger 
values of d while generally decreasing the overlap along 
the way. For the overlap with the Halperin 331 wave- 
function it is clear that there is no region in parameter 
space that produces a sizable overlap in the S'AS'-basis. 
For the layer-basis 331, the overlap is large for moderate 
layer separation (d ~ 1 — 2) and Ag^g. The position 
of the maximum overlap is relatively constant as a func- 
tion of Ap until Ap ~ 0.1 when the maximum overlap 
shifts to higher values of d and decreases markedly in 
magnitude. The expectation values of the pseudo-spin 
operators (Fig. 11:) mirrors the results of the overlap 
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calculations. 

For the FQHE gap shown in Fig. |TT| d we see the famil- 
iar picture found in Ref. El We hrst reiterate that result. 
For Ap = 0, the FQHE gap is largest along a "ridge" 
in AsAS-d space and the overlap between the S'AS'-basis 
Halperin 331 state is largest in the region of phase space 
where the system is two-component and the overlap be- 
tween the SAS-basis Pf state is largest in the region of 
phase space where the system is largely one-component. 
The FQHE gap "ridge" seems to function as the phase 
boundary between the two states in the quantum phase 
diagram, however, the maximum FQHE gap is slightly 
on the Halperin 331 side of the phase diagram and the 
experimentally observecP^HM FQHE at v = 1/2 is of that 
nature. This is strongly thought to be the case because of 
recent work by Storni et a/!^that showed that in the low- 
est Landau level the FQHE gap for one-component sys- 
tems vanishes in the thermodynamic limit-even though 
it is non-zero in our finite size calculation. Thus, our 
non-zero gap in the Pf region of the phase diagram for 
the lowest Landau level is most likely a finite size effect. 

For non-zero Ap the above picture changes. The 
"ridge" is weakened and a maximum in the FQHE gap 
starts to appear for large A$as an d large d until even- 
tually weakening further. In the large Ap limit we see a 
SAS'-basis Pf state for large Ag^g and non-zero d and a 
layer-basis Pf for moderate Asyis and very small d. How- 
ever, we caution that the FQHE gap is globally weakened 
upon inclusion of Ap and taking into account recent re- 
sults^, an experimental system would most likely not 
exhibit the FQHE in that region of parameter space. 

To summarize our results from the calculations on the 
sphere in the lowest Landau level, we find a robust layer- 
basis 331 state that has high overlap with the exact 
ground state and dominant gap in the phase diagram. 
Furthermore, we find a one-component state that has 
some properties of the Moore-Read Pfaffian, but it is 
likely to show up as a compressible state in experiments. 
These conclusions are in agreement with the results of 
Ref. 83 (and consistent with Ref. 58) where transi- 
tions between 331 state, Pfaffian and Composite Fermion 
Fermi sea were studied in a bilayer model with tunneling 
Asas (in the layer-basis) using exact diagonalization and 
effective mean field BCS theory of Read-GreerPll. There 
it was found that the increase of inter-layer tunneling 
converts the 331 state into a Composite Fermion Fermi 
sea because the effective chemical potential of "even"- 
channel electrons becomes very large. Similar analysis is 
expected to apply when Ap is large and the system is 
viewed in S'AS'-basis. 



B. Wide-quantum- well 

We now turn to the wide-quantum-well (WQW) model 
which largely mirrors the results presented for the bilayer 
model. However, there are some differences between the 
two models which we point out below. 



Fig. [12] shows the overlap between the Moore-Read 
Pfaffian wavefunction (in both the S'AS'-basis and the 
layer-basis) and the exact ground state of the WQW 
model as a function of Ag^g and Ap for a few values 
of the WQW width W . This figure should be compared 
with Fig. [6] for the bilayer model. Qualitatively, the two 
models produce very similar results. Of course, a layer 
separation of d in the bilayer model is not equivalent to 
the WQW width W and any similarities between the two 
at d ~ W is coincidental. However, the behavior for the 
bilayer model as d increases is the same as the behav- 
ior for the WQW model as W increases. That is, for 
small W — 0.4, the overlap between the exact ground 
state and the one-component Pf in the layer-basis be- 
comes large when Ap is increased since the system is be- 
ing driven to be more one-component in the layer sense. 
When As as is increased, the overlap between the exact 
ground state and the Pf in the S'AS'-basis becomes large 
since the system is being driven to be one-component in 
the SAS-basis. 



Note, however, that for W = 0.4 (Fig. 12 1) the sys- 
tem is not as symmetric upon As as Ap as it was 
for the d — bilayer model case where the model was 
actually SU(2) symmetric. This is because the WQW 
model is not SU(2) symmetric for small W because even 
as W becomes very small the Coulomb energy between 
electrons in the S and AS states is never equal, i.e., the 
AS state always has a node and therefore higher kinetic 
energy. That being said, the models produce very similar 
results. 

In Fig. [13] we show the overlap between the exact 
ground state of the WQW model and the Halperin 331 
state. This time, as opposed to Fig. [7] we only show the 
overlap with the 331 state written in the layer-basis since, 
as we learned previously by studying the bilayer model 
results, and by confirming this with the WQW model, the 
overlap between the ground state of the WQW model and 
the S^S-basis 331 state is always nearly zero, hence, we 
do not bother to show these results explicitly. For the 
layer-basis 331, however, we again see qualitatively sim- 
ilar behavior compared to the bilayer model. For small 
W the overlap with the 331 state is very small and it can 
be increased by increasing W. One interesting difference 
between the results of the two models is that for W — 5 
the maximum overlap with 331 has a maximum for non- 
zero Ag^g for the WQW model. Furthermore, non-zero 
charge imbalance Ap eventually destroys the Halperin 
331 state by driving the system to be one-component, 
however, the 331 state is remarkably robust to charge 
imbalancing. 

Similar to the bilayer model, the pseudo-spin expec- 
tation value for the WQW model in Fig. [14] mirrors the 
behavior of the overlaps. When the system is largely 
one-component, either in the layer- or 5A5-basis (large 
value of (N R - N L )/2 or (N s - N A s)/2, respectively), 
the overlap with the appropriate basis Pf state is large. 
When the system is two-component the overlap with the 
331 is large. In investigating the pseudo-spin expecta- 
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(a) lll, nfco.4 



(b) LLL, W=2.0 




(b) LLL, W=2.0 



(c) LLL, Hfc4.0 
0.3 



FIG. 12. (Color online) Lowest Landau level: Wavefunction 
overlap between the Moore-Read Pfaffian wavefunction in the 
layer-basis (top panel) and the S'AS-basis (lower panel) and 
the exact ground state for the WQW for (a) W = 0.4, (b) 
W = 2, (c) W = 4 and (d) W = 5. 



(a) LLL, W=0.4 
0.3 



(b) LLL, W=2.0 
0.3. 




(c) LLL, W=4.0 
0. 



FIG. 13. (Color online) Lowest Landau level: Wavefunction 
overlap between the Halperin 331 wavefunction in the layer- 
basis and the exact ground state for the WQW for (a) W = 
0.4, (b) W = 2, (c) W = 4 and (d) W = 5. 



tion value for the WQW model we see the most serious 
discrepancy between the two models-this effect was also 
seen previously in Ref. |31 For W > 2, and relatively 
small values of Ag^s, the ground has a negative value 
of (Ns — Nas)/%- This means that the electrons occupy 
the AS state compared to the S state. If we recall the 
overlap between the exact ground state and the Moore- 
Read Pfaffian in the .S'AS'-basis (lower panel of Fig. 12 ) 




FIG. 14. (Color online) Lowest Landau level: (Pseudo-spin) 
expectation value of the exact ground state (for the WQW 
model) of (N R -N L )/2 (top panel) and (N s - N A s)/2 (lower 
panel) for (a) W = 0.4, (b) W = 2, (c) W = 4 and (d) W = 5. 
Note that the WQW model always breaks SU(2) symmetry 
even for small W. 
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FIG. 15. (Color online) Lowest Landau level: FQHE energy 
gap for the WQW model as a function of W and Asas with 
(c) Ap = (left column) and (d) Ap = 0.1 (right column). 
Also shown is the overlap between the S'ylS'-basis Moore-Read 
Pfaffian (top panel) and the layer-basis Halperin 331 (bottom 
panel) wavefunctions and the exact ground state for (a) Ap = 
and (b) Ap = 0.1. 
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FIG. 16. (Color online) Ratio of the few strongest, inter- 
layer and intra-layer (bilayer) pseudopotentials V£ /V£? 
derived using the WQW modeP as a function of width W. 
Note the saturation of the ratios for large widths W, illus- 
trating the limit of validity of the WQW to describe bilayers. 



the overlap is large for large As^s and this behavior mir- 
rors the large value of (N$ — Nas)/% shown in the lower 
panel of Fig. 14 However, the one-component Moore- 
Read Pfaffian we consider pairs electrons in the S state 
(shown schematically in Fig. [4]o) . We have not checked 



the overlap with the one-component Pf written such that 
the pairing occurs between electrons in the AS state, 
thus, adding another version of the Pf. By examining 
the value of (N s - N AS )/2 for W > 2 in the small A SAS 
region of phase space we would expect a reasonable value 
of the overlap with a one-component Pf pairing electrons 
in the AS state. 

The WQW model is more general than the usual bi- 
layer model and the latter can be derived from it (see Sec- 
tion IV-B of Ref. HJfor a detailed discussion). However, if 
the well width W becomes very large, it is unjustified to 
restrict the model to only the lowest two subbands (and 
the inclusion of higher subbands cannot be treated ex- 
actly because of the computational complexity) . On the 
other hand, as we map the WQW Hamiltonian to an ef- 
fective bilayer, we find that the effective bilayer distance 
d saturates for large W. To see this, we can calculate the 
Haldane pseudopotentials V^ lff2<73Cr4 , where m is the rela- 
tive angular moment of the two electrons and a — S, AS. 
If we denote by J^ a ^ a 'i a ^i(^ ^ ne effective interaction in 
the plane, we have 

Fai(T 2 a z a i {T) = / dzx dz 2 (r\a 1 )(z 1 )(r\a 2 )(z 2 ) 
Jo Jo 

<r|(T3>(*i)<r|(7 4 >(*3) . (8) 
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where C m is the Laguerre polynomial and the Fourier 
transform j^ 1<T2 °' 3<T4 can be evaluated analytically. Fur- 
thermore, we can construct linear combinations^ of 
j/CTiCT 2 o- 3 o- 4 ^0 get the effective bilayer pseudopotentials, 
V™ tra and V™ tCT . We plot the ratios of the few strongest 
T4 ntra and V™ tcI in Fig. Il6] as a function of W. Notice 



The Haldane pseudopotentials for the effective interac- 
tion, written on the disk for simplicity (for the lowest 



that the limits saturate for large W , indicating that the 
inter-layer repulsion decreases very slowly with respect 
to intra-repulsion for larger W, thus suggesting that the 
model becomes unrealistic in this regime. Also that the 
^pseudopotentials involving a node in the z-wave function 
(the ^45 single particle wavefunciton has one node) can 
be smaller in absolute value than those without a node, 
indicating a possibility for some W to have a depopu- 
lation of the lowest subband and negative polarization 
(N s — N AS ), a s seen in the data. 

Fig. [15] shows the FQHE energy gap for the WQW 
model for Ap = and Ap = 0.1, respectively. For Ap = 
0, we see similar behavior to the results of the bilayer 
model and shown previously^, i.e., Fig. Ill for Ap = 
0. Of course, there are quantitative differences between 
the two models but the FQHE energy gap still shows a 
prominent "ridge" as a function of As^g and W. The 
ridge marks the transition line between the 331 and the 
Pfaffian, or perhaps compressible states. Also, note that 
the maximum FQHE gap (the "ridge" ) is slightly on the 
Halperin 331 side of the quantum phase diagram if the 
phase boundary is taken to be the place at which the 
331 overlap is larger than the Pfaffian overlap, i.e., the 
331 phase is the region where the Halperin 331 overlap 
is larger than the Pfaffian and vice versa, cf. Figs. [T5^ 
and b. When the charge imbalancing term is increased 
to Ap = 0.1 the FQHE gap is markedly reduced, as it 
was in the bilayer model. Again, as for Ap = the 
maximum gap is in the Halperin 331 part of the phase 
diagram. Interestingly, there appear to be two "ridges" 
forming for the Ap — 0.1 situation in the FQHE gap 
and both "ridge" maxima are mirrored peaks in the 331 
overlap-although, note that the overlaps (both 331 and 
Pf) for the Ap = 0.1 situation are never very large and in 
a real experimental system the phase would most likely 
have been taken over by some other competing phase by 
the time the charge imbalancing strength reaches Ap = 
0.1, namely a striped phase or, perhaps, a (Composite 
Fermion) Fermi sea. 



1. Torus geometry 

In this section we consider the wide-quantum-well 
model in the torus geometry^HSSI w hich has one distinct 
advantage to the spherical geometry (of course, at the 
cost of other disadvantages). Recall in the spherical ge- 
ometry, the filling factor is defined as Urrijv-)-oo N/N^ and 
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FIG. 17. (Color online) Lowest Landau level: Energy spectra (given as the energy AE with respect to the ground state) 
for different pseudo-momenta (the different levels) with (a) Asas = Ap = as a function of WQW width W, (b) W = 6.5 
and Ap = as a function of Asas and (c) W = 6.5 and Asas = when the imbalance Ap is increased from zero. The 
pseudo-momenta corresponding to the Halperin 331 state and Moore-Read Pfaffian state ((K x ,K y ) = (0,0) (solid red circle), 
(0,4) (solid blue square), (4,0) (solid green diamond) and (4,4) (orange star) and (K x ,K y ) — (0,4), (4,0) and (4,4), respectively. 
The black dashes correspond to states with pseudo-momenta not belonging to either the 331 or Pf states. 



the relationship between TV and N$ for the Moore-Read 
Pfaffian and Halperin 331 is Nj, — 2N — 3 where the "-3" 
is a consequence of the curvature of a spherical surface. 
Other competing phases for a half-filled Landau level, 
such as a (Composite Fermion) Fermi sea, have a different 
shift and comparing different states on an equal footing 
requires great care and extrapolation to the thermody- 
namic limit. However, on the torus, the filling factor is 
uniquely defined as simply N/N^ and, therefore, different 
states can be directly compared for finite size systems-of 
course, this does not mean a conclusion made for a finite 
sized system will maintain as the thermodynamic limit is 
approached. (We note that it makes no sense to consider 
overlaps extrapolated to the thermodynamic limit since 
they trivially vanish.) 

The torus geometry is generally defined using a do- 
main with sides of length a and b with a ^ b. The aspect 
ratio of the toroidal system is t = a/b. The magnetic 
field does not allow the use of the usual translation op- 
erators but many-body states can be founcP^ using the 
so-called magnetic translation operators with conserved 
pseudo- momenta (K x , Ky). The pseudo- moment a belong 
to a Brillouin zone with (K x = 2irs/a,K y = 2nt/b) with 



s, t = 0, . . . , Nq — 1 with N being the greatest common 
divider of N and N$ . 

Different FQHE states on the torus can be specified by 
their ground state (topological) degeneracy. Generically, 
for a Landau level filling factor of p/q there is always a 
center-of-mass degeneracy equal to q which is invariant 
to the form of the Hamiltonian and, thus, of no physical 
significance- we shall ignore it. Also there can be addi- 
tional degeneracies that occur at special points in the 
Brillouin zone, such as at certain points in a hexagonal 
Brillouin zone, and we will consider these trivial. Fi- 
nally, there can be degeneracies that are related to the 
specific topological nature of certain ground states and, 
therefore, non-trivial. For the two-component (Abelian) 
Halperin 331 stateP we expect a quadruplet of states 
(up to the center of mass degeneracy, which in our case 
is 2) one of which belongs in the (K x , K y ) = (0, 0) sector 
and the remaining three are at the Brillouin zone corners 
(K x ,K y ) = (0,N /2), (N o /2,0), and (N /2,N /2). The 
non- Abelian Moore-Read PfaffiarJ^^I nas only a three- 
fold degeneracy^ of (K x ,K y ) = (0,N /2), (N /2,0), 
and [No/2, N /2), i.e., the (K x ,K y ) = (0,0) state is 
missing. Compressible states, such as the (Composite 
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Fermion) Fermi sea, generally do not have clearly defined 
degeneracies-they may have accidental degeneracies that 
are strong functions of the aspect ratio r, particle num- 
ber N, or other Hamiltonian parameters. The exact de- 
generacies are expectations based on the (analytic) form 
of the variational ansatz wa vefunc tions and their respec- 
tive conformal field theories^^2l. The ground state(s) 
of an actual H amiltonian will not display exact degen- 
eracies^EnHoHni (perhaps they do in the thermodynamic 
limit?) but the ground state(s) should qualitatively show 
the ground state topological degeneracy corresponding to 
the variational ansatz if they are to be thought of as be- 
ing in that "phase" . Many states are sensitive to changes 
in r while others are not. One should really investigate 
the properties of the system with regard to changes in r, 
however, in the present case, a somewhat involved ana- 
lytical calculation regarding the "background" charge is 
needed, and so we focus here on a fixed aspect ration of 
t = 0.97. 

All results on the torus correspond to a system of 
N = 8 electrons in the half-filled lowest Landau level- 
v = 1/2. Fig. 17i shows the energy spectrum of the 



low-lying states for the WQW model in the absence of 
any tunneling terms, i.e., Ag^g = Ap — 0, as a function 
of the WQW width W. For each W, the lowest-lying en- 
ergies are plotted relative to the ground state and shape- 
and color-coded to emphasize which states belong to 
which pseudo-momenta sectors expected when consider- 
ing the two-component Halperin 331 and one-component 
Moore-Read Pfaffian wavefunctions. Namely, we track 
the pseudo-momenta sectors (K x , K y ) = (0, 4), (4, 0) and 
(4,4) for the Pf with the addition of (K x ,K y ) = (0,0) 
for the 331 state. For large widths W > 4.5, the spec- 
trum is characterized by a large energy gap separat- 
ing a high-energy (quasi-)continuum of states from a 
low-energy, nearly degenerate manifold of states at the 
pseudo-momenta corresponding to the 331 state. As the 
width is decreased towards zero the (K x ,K y ) — (0,0) 
state that specifies the Halperin 331 state from the 
Moore-Read Pfaffian state goes up in energy joining the 
continuum. However, the (K x ,K y ) = (0,4) and (4,0) 
states of the Moore- Read Pfaffian also rise in energy into 
the continuum while at the same time a few states from 
the continuum drop down in energy mixing with states 
belonging to the 331 or Pf states. Finally, at very small 
W < 1 there is a single ground state of (K x , K y ) = (4, 4) 
and a large energy gap. However, there is no three- 
fold degeneracy characteristic of the Moore- Read Pfaffian 
state present. 

In Fig. [T7p we fix W = 6.5 and Ap — and vary the 
inter-layer tunneling strength Ag^g. For small A^g = 
we are clearly in the Halperin 331 part of the phase di- 
agram since the spectra has a quasi-four-fold degeneracy 
characteristic of the 331 phase separated from the contin- 
uum by a large energy gap. Upon increasing As^s in an 
attempt to drive the system into the Moore-Read Pfaf- 
fian phase we see that while the (K x , K y ) = (0, 0) state of 
the 331 phase rise in energy and joins the continuum, the 



(K x , K y ) = (4, 4) state of the Pfaffian phase goes with it 
and a state from the high-energy continuum drops down 
into a quasi-three-fold degeneracy. However, this quasi- 
three-fold degeneracy does not contain the right pseudo- 
momenta to describe the non-Abelian Moore-Read Pfaf- 
fian phase. 

Lastly, in Fig. [l7p we again set W — 6.5 but now fix 
Asas = and vary the charge imbalance Ap. As before, 
for small Ap we see a clear signature of the Halperin 
331 phase. However, as Ap is increased the energy gap 
essentially collapses giving way to a phase that would 
most likely not exhibit the FQHE. 

These conclusions corroborate our previous work us- 
ing the spherical geometry and agree with the previ- 
ous studySS of bilayer model on the torus for v = 1/2. 
Namely, the Halperin 331 state is a good ansatz for the 
FQHE in the right parameter regions for bilayer and 
WQW systems. However, when the system is driven to 
be one-component in the hopes of producing a FQHE 
described by the Moore-Read Pfaffian state the Hamil- 
tonian details and Haldane pseudopotentiaP^ values are 
such that the Moore-Read Pfaffian phase loses out to a 
different, most likely, non-FQHE state such as a striped 
phase or (Composite Fermion) Fermi sea. 



V. RESULTS: SECOND LANDAU LEVEL 

In this section we present results of our calculations 
for the second Landau level, that is, bilayer FQHE at 
v = 5/2. (Worlds has been done considering bilayer 
FQHE where the total filling factor is v = 5/2 + 5/2 = 5 
but it is unrelated to our work.) We are neglecting Lan- 
dau level mixing and considering the electrons occupy- 
ing the second Landau level to be spin-polarized. Oper- 
ationally, we have projected the half-filled electrons in 
the second Landau level into the lowest Landau level 
using the Haldane pseudopotentials^. These effective 
pseudopotentials take into account the screening of the 
Coulomb interaction that occurs between the electrons 
in the SLL due to the (taken to be) inert electrons in 
the LLL. We will not detail the procedure of using Hal- 
dane pseudopotentials in the FQHE as this procedure has 
been given in many places^El] Also, we (Peterson and 
Das SarmeP) have recently studied the FQHE in the SLL 
without the presence of a charge imbalancing term and 
will compare our results here extensively with the ones 
given previously. 

Before we tackle the results for the second Landau level 
we briefly note that bilayer systems in higher Landau 
levels are quite subtle and non-trivial and actually pro- 
vide some conceptual difficulty. As recently discussed 
in Ref. 5 ( it is not obvious what happens to the elec- 
trons in the lowest spin-up and spin-down Landau levels 
when driving a one-component system at v = 5/2 to a 
two-component (bilayer) system at total filling v = 5/2. 
For this work, however, we take the conceptually well- 
defined and straightforward solutiorP: we assume full 
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FIG. 18. (Color online) Second Landau level: (a) Wavefunction overlap between the exact ground state of the bilayer Hamil- 
tonian and the Moore-Read Pfaffian written in the layer-basis (left column) and the S^S-basis (right column) shown as a 
function of inter-layer tunneling Asas and charge imbalance Ap for different values of layer separation d and zero individual 
layer thickness w = 0. (b) Same as (a) but for the Halperin 331 wavefunction. (c) Pseudo-spin expectation value or, more 
physically, the expectation value of [Nr — Nl)/2 (left column) and (N$ — Nas)/2 (right column) as a function of Ag^s and 
Ap. (d) The FQHE energy gap for the exact bilayer Hamiltonian. 



spin-polarization, hence the system is essentially spin- 
less with each Landau level having only one spin-index. 
Then the v = 5/2 (balanced) two-component system is 
equivalent to one where each layer has 1 + 1/4 filling 
with the lowest spin Landau level being completely full 
and the second Landau level being 1/4 full. In this way, 
the incompressible FQHE states (Halperin 331 or Moore- 
Read Pfaffian) form completely in the second Landau 



level. This "solution" provides a completely well-defined 
mathematical problem. Of course, the real physical sys- 
tem, i.e., bilayer FQHE systems in higher Landau levels, 
could be more complicated and produce rich physics. In 
fact, theoretically, the full solution is out of reach for any 
conceivable computer-the Hilbert space is too vast when 
including many (or at least three) Landau levels, spin-up 
and spin-down, and layer index. Thus, theoretical and 
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(a) SLL, d=S> (b) SLL, rf=l 




FIG. 19. (Color online) Second Landau level: Wavefunction 
overlap between the Moore-Read Pfaffian wavefunction in the 
layer-basis (top panel) and the S'AS-basis (lower panel) and 
the exact ground state for (a) d — 0, (b) d = 1, (c) d = 2 and 
(d) d = 4. 





FIG. 21. (Color online) Second Landau level: (Pseudo-spin) 
expectation value of the exact ground state of (Nr — N^)/2 
(top panel) and (Ns — Nas)/2 (lower panel) for (a) d = 0, 
(b) d = 1, (c) d = 2 and (d) d = 4. 



(a) SLL, d=0 (b) SLL, ifcl 




FIG. 22. (Color online) Second Landau level: FQHE energy 
gap for (a) d = 0, (b) d = 1, (c) d = 2 and (d) d = 4. 

experimental efforts in studying bilaycr FQHE systems 
in higher Landau levels is likely a fertile ground for new 
discoveries. 



FIG. 20. (Color online) Second Landau level: Wavefunction 
overlap between the Halperin 331 wavefunction in the layer- 
basis (top panel) and the SAS'-basis (lower panel) and the 
exact ground state for (a) d = 0, (b) d = 1, (c) d = 2 and (d) 
d = 4. 



A. Bilayer 

In Fig. [18] we show overlaps between the exact ground 
state and the layer-basis and S'AS'-basis Pf and 331 wave- 
functions, pseudo-spin expectation values, and the FQHE 
gap for the second Landau level. The difference between 
the LLL and the SLL are subtle but important. 
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First we focus on the overlaps between the Pf and 331 
states and the exact ground state in both the layer- and 
S'AS-bases, cf. Fig. [l8^-b. The main difference between 
the results in th e LLL and the SLL is that, as has been 
shown previousl y ^ 5 * 59 * 60 ! for the case of zero charge imbal- 
ance Ap = 0, the overlap between the exact state and the 
Pf state(s) is higher in the SLL than it is in the LLL. In 
the SLL, the Pfafhan overlap is approximately ~ 0.97 at 
its largest while in the LLL it is approximately ~ 0.9 at 
its largest. This can be seen more clearly in Fig.|l9|where, 
compared to the results of the LLL, the overlap with the 
Pf state is higher in the SLL. (Again, for d = the system 
is SU(2) symmetric and the layer-basis and S'AS'-basis 
results are related via a pseudo-spin rotation.) We note 
that by increasing the width of the individual quantum 
well w, the Mo ore- Re ad Pfaffian overlap can be increased 
to nearly unity^2Ml This is a consequence of the differ- 
ences between the Haldane pseudopotentials correspond- 
ing to the LLL versus those of the SLL with the SLL 
pseudopotentials being more amenable to pairing into a 
p-wave paired BCS stat e (of Composit e Fermions), i.e., 
the Moore-Read Pfaffia ri 42 ' 50 ' 59 ™ 9 ^. 

There is another striking difference between the results 
in the LLL versus the SLL when considering the overlap 
between the Halperin 331 state. In the LLL, the 331 
state provides a very good description of the FQHE for 
two-component systems when the separation d (or wide- 
well width W for the wide-quantum-well) and the tun- 
neling strengths are appropriate-this is evident by the 
value high value of the overlap that is obtained (~ 0.99). 
In the SLL, the qualitative behavior of the overlap as the 
system parameters are varied is similar to the LLL but 
the overlap does not obtain as high a value (only ~ 0.8 
for d ~ 1 and Ag^g = Ap = 0). This, again, can be seen 
more clearly in Fig. [20] 

Looking at the pseudo-spin expectation value (S z ) = 
{N s - N AS )/2 and (S z ) = (N R - N L )/2 in Fig. g (as 
well as Fig. 21 ) we see very little difference between the 
lowest and second Landau levels. In fact, visually it is 
difficult distinguishing the two. 

Lastly, we consider the FQHE gap (Fig 



18 i and 



Fig. 22 1, which does display some behavior that is quali- 



tatively different from that of the LLL. This qualitatively 
different behavior manifests itself for values of layer sep- 
aration d > 2. In this region of parameter space we see 
that for increasing values of Ap the FQHE gap has a max- 
imum for a wide range of Asyis. This is in stark contrast 
to the FQHE gap in the LLL (cf. Fig. [5]i) where the gap 
shows some indication of increasing for increasing Ap but 
not to a maximum. 

Fig. 23 shows the overlaps, pseudo-spin expectation 



values, and FQHE gap as a function of layer separa- 
tion d and tunneling strength A 5^45 for several values of 
charge imbalance tunneling strength Ap. Qualitatively, 
the overlaps (Fig. 23 1 and b) are similar in the SLL and 
LLL. However, comparing them side-by-side, the overlap 
with the Halperin 331 state achieves a higher value in 
the LLL than it does in the SLL. The opposite is true for 



the overlap with the Moore-Read Pfaffian state where the 
overlap achieves a higher value in the SLL than it does 
in the LLL. These differences are manifest even while 
the pseudo-spin expectation values (Fig. [23];) in the SLL, 
compared to the LLL, cannot visually be distinguished. 
Thus, while the values of the Haldane pseudopotentials 
in the SLL are not different enough from the LLL to 
change pseudo-spin expectation values, i.e., the system 
is one-component or two-component at almost exactly 
the same place in parameter space in both the LLL and 
SLL, the differences arc enough to change the character 
of the overlaps and the FQHE energy gap. 

Fig.[23p! shows the FQHE energy gap, and again, there 
is a marked difference between the result in the LLL ver- 
sus SLL. Similar to the higher overlap value between the 
exact ground state and the Moore- Read Pfaffian state in 
the SLL, the FQHE energy gap has a maximum in the Pf 
region of the approximate quantum phase diagram espe- 
cially along the As as axis for small d. (The opposite is 
true in the LLL-the FQHE energy gap is largest in the 
Halperin 331 region of the quantum phase diagram.) We 
see this behavior also in the wide-quantum-wcll results. 



B. Wide-quantum- well 

In Fig. [24] we plot the overlap between the exact ground 
state of the WQW model and the Moore-Read Pfaffian 
wavefunction written in the layer-basis (top panel) and 
5A5-basis (lower panel). This figure is, of course, simi- 
lar to, and should be compared to, Fig. [T9] showing the 
same thing for the bilayer model. Qualitatively, the same 
behavior is manifest in the two models. As the WQW 
width W is increased it takes a larger and larger tunnel- 
ing strength to increase the overlap with the Pf. 

The overlap between the ground state and the Halperin 
331 state is shown in Fig. [25] Again, we only show the 
overlap with the Halperin 331 state written in the layer- 
basis since the overlap with the SAS-basis 331 state is 
almost zero throughout parameter space. The results 
are qualitatively similar to those of the bilayer model. 
However, for W — 2 (Fig. 25 3) there is a slight maxi- 



mum in the overlap as Asas is increased from zero at 
Asas ~ 0.03. This is an interesting result that is harder 
to understand-increasing Asas drives the system to be 
one-component and one would surmise that increasing 
this parameter would lower the overlap with the 331 state 
monotonically. For W — 4 we see that two peaks in 
the overlap appear; one for small values of Ag^g and 
another for a larger value with a strong minimum in- 
between; behavior is particular to the WQW model. Sim- 
ilar to the lowest Landau level, the expectation value of 
(Ng — Nas)/2 for W > 2 is negative for small Ag^g and 
Ap. In fact, this negative value means that the system 
prefers to have more electrons in the AS state than the S 
state. However, the system is still more one-component 
in those regions, such as the region in Fig. [26p along the 
Ap = axis near Asyts ~ 0.3, which correspond to a 
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FIG. 23. (Color online) Second Landau level: Same as Fig. [18] except that all plots are shown as a function of inter-layer 
tunneling Asas and layer separation d for a number of different charge imbalances Ap and zero individual layer thickness 
w = 0. 



maximum in the overlap with the 331 state. At any rate, 
this behavior marks a qualitative difference between the 
bilayer and WQW models. 

Lastly, we show the FQHE energy gap for Ap = 



in Fig. 27 3 along with the overlap between the SAS- 
basis Moore-Read Pfaffian state (top panel of Fig. [27^ ) 
and the layer-basis Halperin 331 state (bottom panel of 
Fig. 27 1), respectively. When the FQHE energy gap is 
non-zero, and the system therefore would be expected to 
exhibit the FQHE, the overlap with the Pfaffian state 
is large. On the other hand, in the region of parameter 
space where the FQHE energy gap is small, the Halperin 



331 state has a larger overlap. Thus, independent of 
our model of choice, the details of the electron-electron 
interaction in the second Landau level, compared to the 
lowest Landau level, make a profound difference. In the 
second Landau level the Pfaffian wavefunction is a good 
description of the physics when the system is largely one- 
component. In the lowest Landau level, the Halperin 331 
is a good description of the physics when the system is 
largely two-component. 

Note that there is a peak in the FQHE energy gap for 
W ~ 2 and Ag^g sa which is well within the Halperin 
331 part of the approximate quantum phase diagram. 
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FIG. 24. (Color online) Second Landau level: Wavefunction 
overlap between the Moore-Read Pfaffian wavefunction in the 
layer-basis (top panel) and the S'AS-basis (lower panel) and 
the exact ground state for the WQW for (a) W = 0.2, (b) 
W = 2, (c) W = 4 and (d) W = 8 as a function of A S as and 
Ap. 
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FIG. 25. (Color online) Second Landau level: Wavefunction 
overlap between the Halperin 331 wavefunction in the layer- 
basis and the exact ground state for the WQW for (a) W = 
0.2, (b) W = 2, (c) W = 4 and (d) W = 8 as a function 
of Asas and Ap. Note that the overlap with the layer-basis 
Halperin 331 state is not shown since it is always nearly zero. 



Even though the overlap with the 331 state is small in 
that region of parameter space it could correspond to a 
two-component Abclian FQHE for bilayer systems in the 
second Landau level at total v = 5/2 that is in the same 
universality class as the Halperin 331 state. 
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FIG. 26. (Color online) Second Landau level: Pseudo-spin 
expectation value of the exact ground state (for the WQW 
model) of (N R - N L )/2 (top panel) and (7V S - N A s)/2 (lower 
panel) for (a) W = 0.2, (b) W = 2, (c) W = 4 and (d) W = 8. 
Note that the WQW model always breaks SU(2) symmetry 
even for small W and that for W > 2 there are negative values 
of (Ns — Nas)/2 meaning the system prefers, in some regions 
of parameter space, to have more electrons in the AS state 
compared to the S state even for positive and non-zero Asas- 
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FIG. 27. (Color online) Second Landau level: (b) FQHE en- 
ergy gap for the WQW model as a function of W and Asas 
with Ap — 0. Also shown is the (a) overlap between the 
SAS-basis Moore-Read Pfaffian (top panel) and the layer- 
basis Halperin 331 (bottom panel) wavefunctions and the ex- 
act ground state for Ap = 0. 
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VI. CONCLUSIONS 

In conclusion, we have investigated the FQHE in two- 
component systems for the half-filled lowest and second 
Landau levels (v — 1/2 and v = 5/2, respectively) as a 
function of both tunneling strengths; inter-layer tunnel- 
ing and charge imbalancing. This work was motivated by 
the recent experimental systems investigated in Refs. 6 
and [7] 

Our main results are as follows: 

(i) The FQHE at v = 1/2 is described by the (Abelian) 
Halperin 331 two-component state which is remarkably 
robust to charge imbalance tunneling. When the sys- 
tem is driven into the one-component region of param- 
eter space, we find that the (non- Abelian) Moore-Read 
Pfaffian state is most likely beaten out by other com- 
peting non-FQHE phases (cf. striped phase or (Com- 
posite Fermion) Fermi sea). The reason we suspect this 
is because even though the overlap between the exact 
ground state and the Pf (in the SAS- or layer-basis, de- 
pending on the region of parameter space, i.e., whether 
Asas is large and Ap is small or vice-versa) is large, 
the FQHE energy gap has a peak that is (slightly) on 
the Halperin 331 side of the approximate quantum phase 
diagram, where the quantum phase diagram determina- 
tion is described in Sec. |IV| This result, along with re- 
cent numerical calculations by Storni et alW* (as well 
as Refs. IFTl and I53"j) . leads us to conclude that the one- 
component v = 1/2 is likely to be a non-FQHE striped 
phase or (Composite Fermion) Fermi sea. In the effective 
BCS description, the increase of tunneling (in the layer- 
basis) or charge imbalance (in the iSAS-basis) leads also 
to the increase of the effective chemical potential of the 
"even" channel which drives the system into a compress- 
ible phased. 

We also find that our calculations are unable to explain 
the recent results of Shabani et alWQ which consider only 
the competing Moore-Read Pfaffian and Halperin 331 
states. In that work, they observe no FQHE at v = 1/2 
when the charge imbalance is zero, an emerging FQHE at 
v = 1/2 for non-zero charge imbalance, and finally, that 
the FQHE is eventually destroyed upon increasing charge 
imbalance. While we find that the Halperin 331 FQHE 
is quite robust to a charge imbalance tunneling term, our 
results would suggest that the FQHE would monotoni- 
cally decrease in strength with increasing charge imbal- 



ance (that is, the measured activation gap would decrease 
in strength, the minimum in R xx would weaken, and/or 
the "quality" of the plateau in R xy would deteriorate)-we 
would expect experiments to observe a FQHE at v = 1/2 
with zero charge imbalance that is eventually destroyed 
upon further imbalancing. The fact that this does not 
happen in the experiment indicates that the observed 
state at large imbalance must be something different from 
the 331 or the Pfaffian state, at least within our model 
calculations. 

(ii) The FQHE in the second Landau level at v = 5/2 is 
most likely the spin-polarized and one-component (non- 
Abelian) Moore-Read Pfaffian state. This state, similar 
to the Halperin 331 state in the lowest Landau level, is 
remarkably robust to charge imbalancing. We also find 
that in regions of parameter space when the system is 
largely two-component, i.e., for small Ag^g and Ap, the 
system might display a FQHE described by the (Abelian) 
Halperin 331 state. This suggests^ the exciting possi- 
bility of experimentally tuning parameters to drive an 
Abelian FQHE state (331) at v = 5/2 into a non-Abelian 
FQHE state (Pf) at v = 5/2 by way of a quantum phase 
transition. 

(iii) One of our surprising theoretical findings is that 
the Halperin 331 two-component bilayer Abelian paired 
FQHE is extremely robust, and survives, not only sub- 
stantial inter-layer tunneling, but also substantial charge 
imbalance. We suspect that this result is quite general 
and possibly also applies to other Abelian paired states^ 
in the same universality class of the 331 state. 

The physics of the bilayer FQHE at v = 5/2 is ex- 
tremely rich with many unanswered and fascination ques- 
tions (not to mention the fact that the bilayer problem 
in higher Landau levels is conceptually difficult and non- 
trivial (see Sec. |V|) that, in our opinion, are awaiting 
experimental answers. 
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